Magnetism in C or N-doped MgO and ZnO: density-functional study of impurity pairs 
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It is shown that substitution of C or N for O recently proposed as a way to create ferromagnetism in 
otherwise nonmagnetic oxide insulators is curtailed by formation of impurity pairs, and the resultant 
C2 spin=l dimers as well as the isoelectronic Nj^ interact antiferromagneticalUy in p-type MgO. 
For C-doped ZnO, however, we demonstrate using the HSE hybrid functional that a resonance of the 
spin-polarized C2 ppn* states with the host conduction band results in a long-range ferromagnetic 
interaction. Magnetism of open-shell impurity molecules is proposed as a possible route to 
ferromagnetism in oxide spintronic materials. 

PACS numbers: 75.50.Pp, 71.20.-b, 71.70.-d 



Ferromagnetie (FM) semieonductors or insulators that 
are free of transition metal or rare earth species are nowa- 
days perceived as a class of systems displaying (f mag- 
netism [iHll . As a common feature, defects with sp-type 
orbital character carry the magnetic moment in these 
systems. One widely discussed scenario for functionaliz- 
ing oxides as dP dilute magnetic semiconductors involves 
the magnetic polarization of valence states by the sub- 
stitution of oxygen with sp-type impurities of lower va- 
lency. To maximize the Hund exchange at the impu- 
rity sites, doping by the 2p elements C or N appears 
most promising. For example, local-density approxima- 
tion (LDA) calculations predicted that group-II oxides 
MO (M=Mg^ Ca, Sr, Ba) become FM when doped with 
C or N 0, m. Moreover, experiments and computa- 
tions demonstrated that C-doped ZnO is FM even above 
room temperature 04^- Typically samples with very 
high concentration of impurities are used. Under these 
conditions, it was believed that the FM interaction be- 
tween the defects is mediated by a partially occupied 
spin-polarized defect band. However, this mechanism re- 
quires considerable fine-tuning: The defect states need 
to be sufficiently extended to mediate the FM coupling 
via Zener's double exchange, yet the band width must be 
small enough for the Stoner criterion to be satisfied [13- 
[lit . Consequently, theoretical results reported so far are 
strongly method-dependent: In LDA or generalized gra- 
dient approximation (GGA) calculations, the 2p defect 
state tends to come out too extended, and hence their 
predictions of ferromagnetism might be overly optimistic. 
On the other hand, calculations using LDA plus Hubbard 
U [1^ or a self-interaction (SIC) correction em- 
phasize Coulomb correlation in the 2p hole states, and 
favor a splitting of the partially occupied impurity lev- 
els. The Coulomb correlation strongly localizes the 2p 
hole states (and the associated magnetic moments) and 
thus, in combination with a local Jahn- Teller distortion. 



impedes ferromagnetism 

In this Letter, wc study impurity pairs that are ex- 
pected to play a significant role at the high impurity con- 
centrations typically used. We demonstrate that C- or 
N-impurities may interact attractively in both MgO and 
ZnO. In contrast to the earlier predictions of a metal- 
lic impurity band, we find that MgOi-a^N^;, after pair- 
ing, becomes an insulator with nonmagnetic or antifcrro- 
magnetically coupled impurity molecules. Next, we ad- 
dress the magnetic properties of the open-shell 'defect 
molecule' C2 in these oxides hosts. Our calculations sug- 
gest that the dP magnetism observed in the C-doped ZnO 
may well be due to a magnetic interaction between these 
species mediated by conduction band electrons. 

We performed LDA calculations for the N- or C-doped 
MgO, using fully relaxed structures with two substitu- 
tional impurity atoms in a 2x2x2 supercell. The full- 
potential augment ed p lane wave plus local orbital code 
WIEN2k was used [20| . The muffin-tin spheres were cho- 
sen to be 2.1 Bohr for Mg and 1.3 for O, N and C, the 
plane- wave cutoff of 15 Ry for the interstitial wave func- 
tions, and a 5x5x5 k-point mesh for integration over the 
Brillouin zone. The Kohn-Sham band gap of bulk MgO 
is calculated to be 5.0 eV by LDA, much smaller than the 
experimental value of 7.8 eV. We note that the impurity 
bands of concern are all located inside the (too small) 
LDA gap (see Figs. [T] and [2]). Thus, the results can be 
interpreted in a physically meaningful way, despite the 
obvious shortcomings of LDA. Note also that our calcu- 
lations using GGA gave almost the same results as LDA. 

We start with the calculations of two substitutional N 
atoms at various oxygen sites in a Mg32 03oN2 supercell, 
for a comparison with the literature. Fig. [TJa) shows the 
density of states (DOS) obtained with the longest N-N 
distance of 7.29 A. It is a FM half-metal (FMHM) as 
reported in the literature This is due to the 2/3 

filled down-spin N 2p impurity band split off from the 
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FIG. 1: (Color online) Total and N 2p DOS for two N atoms 
with the largest (a) or smallest (b) possible separation, in 
a Mg3203oN2 supercell. The Fermi level is set to zero, (c) 
Charge density contour plot for case (b) (0.1-0.8 e/A ). The 
dashed-line square marks the supercell in the ab plane, (d) 
Schematic energy level diagram of a N-N pair. Except for the 
ppa antlbonding state, all levels are fully occupied, consistent 
with (b). 



fully occupied up-spin band via Hund exchange. After 
placing the second N atom at the nearest neighbor site, 
the total energy is lowered by 0.4 eV compared to the well 
separated N impurities [Fig. [UJa)] and the N-N distance 
becomes 0.65 A shorter than the 0-0 distance [Fig.[TJc)]. 
Fig. [TJd) shows a schematic level diagram of the N-N 
pairing interaction. The ppa bonding state and the ppn 
states can accommodate all the ten electrons, and the 
ppa antlbonding state is empty; thus the pairing leads to 
a nonmagnetic insulating state, see also Fig. 1(b). 

Now we turn to the C-doped MgO. For the Mg3203oC2 
supercell with the longest C-C distance, our LDA calcu- 
lation gives the same FMHM solution (due to the 1/3 
filled C-2p down-spin impurity band, not shown here) as 
reported in the literature However, the pairing in- 
teraction for C in MgO is even more significant than for 
N: The C-C distance is 1.47 A, being only half of the 
0-0 distance of 2.97 A, see Fig. Mjo). The correspond- 
ing bonding energy gain is large, being 3.2 eV compared 
to well-separated C impurities. This much stronger C-C 
pairing effect can be understood by noticing that there 
are two electrons less in the ppn antlbonding level of a 
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FIG. 2: (Color online) (a) Total and C 2p DOS for two 
nearest-neighbor C impurities in a Mg3203oC2 supercell. It 
is an insulator with a molecular spin=l. (b) Charge density 
contour plot (0.1-0.8 e/A ). (c) Schematic energy level dia- 
gram of a C-C molecular ion, isoelectronic to an O2 molecule. 
The nearly degenerate ppn antlbonding state is half-filled and 
is in a spin triplet, consistent with the DOS in (a). 



C-C pair than in N-N, see Figs. [2{c) and[TJd). It is not 
surprising that the C-C pair has a molecular spin=l, be- 
cause it is isoelectronic to the free O2 molecule, the spin 
triplet ground state of which is well known. Similarly, 
a molecular spin was found in N^" [2l| and O^ 0, 



molecular ions being part of ionic compounds. As seen 
in Fig. HJa), the C-C dimer has a huge ppa bonding- 
antibonding splitting of more than 6 eV, as well as a 
sizable ppn bonding-antibonding splitting of 3 eV. It is 
important to note that the formal ppn doublet has a crys- 
tal field splitting of 0.4 eV due to atomic relaxations, but 
this is still smaller than the exchange splitting of 0.8 eV. 
As a result, the C2 molecular spin=l is preserved in the 
calculations, which show that indeed the spin triplet is 
by 0.19 eV more stable than a spin singlet. Moreover, 
there is a small gap of about 0.4 eV between the spin-up 
and down ppn antlbonding states. Thus, an antiferro- 
magnetic (AF) superexchange is expected between those 
spinal C2 molecules, and our calculations show that 
the AF coupling between two most distant C2 molecules 
in a Mg32 03oC2 supercell is more stable than the FM 
by 20 meV. (We note that this value would be some- 
what reduced if in the calculations the true band gap 
rather than the smaller LDA gap and a larger supercell 
were used). Such an AF coupling was also found in the 
molecular magnets SrN [2l[ and Rb406 [11]. In other 
words, our calculations demonstrate that substitutional- 



3 



C-doped MgO, in the absence of other defects, is expected 
to be an insulator, and has impurities with a molecular 
spin=l and finite AF coupling. Thus, the present re- 
sults are in strong contrast to previous work assuming 
single impurities, predicting a FMHM by LD A m or an 
insulator with a local atomic spin=l by SIC jlTl llSj. 

Next, we address the question whether the presence 
of intrinsic defects or co-doping could possibly change 
the position of the Fermi level Ep and thus induce 
magnetism. Formation energies of single impurities 
and impurity pairs in various charge states were cal- 
culated in the presence of a compensating background 
charge 2J], using a 128- atom supercell for MgO. An 
electrostatic correction for the finite supercell size was 
included [lEl- These GGA calculations were performed 
using the projector-augmented-wave method [26| as im- 
plemented in the VASP code [l^]. The results shown in 
Fig. [3] confirm our negative conclusions regarding mag- 
netism in MgO; Single impurities are energetically favor- 
able only in n-typc MgO, but the resulting and 
species are closed-shell and hence nonmagnetic. In p- 
type MgO (low Ep values) we find that the formation of 
single impurities is energetically very costly, and pairing 
is exothermic. The thermodynamically stable species (cf. 
Fig. [3]) are and , carrying a molecular spin =1 and 
= i, respectively. Notably the Nj"*" becomes more stable 
and magnetic by loosing two antibonding spin-down ppir 
electrons in p-type MgO, and the N-N distance is re- 
duced to 1.29 A. Being isoelectronic to neutral C2 (cf. 
Fig. 2), the insulating band structure of N2'''/MgO (not 
shown here) results in an AF superexchange, as discussed 
above. For C/MgO, the stable Cl+ , Ct and C§ species 
(cf. Fig. 131) carry a molecular spin =0, =^ and =1, re- 
spectively, but again AF coupling (if any) is expected due 
to the samples being band insulators. Thus, the route to 

-magnetism envisaged previously, namely via 2p-holes 
in an impurity band formed from diluted impurities [3~ 
0, [13413, appears not well founded, as isolated impurities 
are far from being thermodynamically stable. The paired 
state, on the other hand, is insulating and hence doesn't 
support FM interactions. 

In view of the above spin=l C2 molecule, we wish 
to find a suitable oxide semiconductor which, after C 
doping, becomes FM. This requires that its band gap 
should be about 3-4 cV, with the bottom of its conduc- 
tion band overlaps the spin-polarized antibonding ppir 
impurity state [see Fig.l^^a)]. Then their level resonance 
[see Fig. 4(a)] will make the bottom of the conduction 
band spin-polarized and result in a long-range FM (28| . 
As such, we propose ZnO as candidate, having a band 
gap of 3.4 eV. First, we calculate the formation energies 
of isolated and paired C in ZnO in a 108-atom super- 
cell with GGA, including an effective Hubbard U=9 eV 
for the Zn 3d electrons to correct the underbinding of 
the 3d states. Again, we find pairing to be preferred (cf. 
Fig. [3]) as the Fermi level runs through the whole band 
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FIG. 3: (Color online) Formation energy of C and N substitu- 
tional impurities at O sites in MgO and ZnO for two isolated 
single impurities (dashed) and one impurity-pair (full lines) in 
different charge states with respect to thermodynamic reser- 
voirs of atomic C and N, and gas-phase O2. The Fermi level 
is measured relative to the valence band maximum. 



gap. In order to explore the consequence of impurity- 
pair levels resonating with the host conduction band, 
we used the screened hybrid functional HSE [sof . 
which improves the description of semiconductors and 
their defect states (3ll - [33| by mixing in a small amount 
(denoted by a) of the Hartrec-Fock exchange to the DFT 
exchange. In improvement over previous LDA and GGA 
studies [3, [III [l3i USE yields, for ZnO, a band gap in 
excellent agreement with the experimental value when 
q:=:0.375 is used [s^l- We have studied the specific case 
of two C-C dimers under the condition of charge neu- 
trality, using the VASP code and a 4x2x2 (Zn32 028 C4) 
supercell together with a 2x4x4 k-point grid. The plane 
wave cut-off energy of 400 eV has been set, and a=0.375 
and the screening factor of 0.2 have been used. After 
relaxation of a nearest-neighbor impurity pair, we obtain 
C-C dimer placed asymmetrically in the ZnO matrix: one 
of the oxygen vacancies is occupied by the dimer, while 
the other vacancy shrinks as the adjacent Zn atoms relax 
inward. The C-C bond length is 1.31 A, even shorter than 
for the C-C dimer in MgO matrix. The two C2 dimers 
are separated by 6.5 A . The C2 molecule in a single va- 
cancy hybridizes considerably with the neighboring Zn 
atoms. Hence, the antibonding molecular states broaden 
into resonances with the conduction band and the spin 
density spreads over the C2 dimers and the neighboring 
Zn atoms, see Fig. 4(b). The spin moment is calculated 
to be 1.4 fiB per C2 dimer (reduced from its formal value 
of 2 (Ib) due to the hybridization and partial delocal- 
ization. It is most important to note that after a full 
relaxation of the supercell having a FM or AF coupling 
between the two C2 dimers, the FM state turns out by 
our HSE calculations to be more stable than the AF state 
by 36 meV, indicating that C-doped ZnO may indeed be 
a FM semiconductor. 

To conclude, we demonstrated by LDA calculations 
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FIG. 4: (Color online) (a) Resonance of the spin-split C2 
impurity levels with the conduction band bottom of ZnO sta- 
bilizes FM order. Solid (dashed) DOS curves standing for the 
switch-on (off) hybridization illustrate the energy lowering, 
(b) Top view into the ab plane of the spin-density isosurface 
of the C doped ZnO in the FM ground state. Big blue (small 
red) spheres stand for Zn (O) atoms. The spin density spreads 
over the C2 dimer and the neighboring Zn atoms. 



that C or N substitutional impurities in p-type MgO 
prefer pairing, which renders N-doped MgO a spin=l 
antiferromagnetic (spin=0 nonmagnetic) insulator with 
the 2+ charged (neutral) N-N pairs, and C-dopcd MgO 
an insulating molecular antiferromagnet. Moreover, we 
find that C-dopcd ZnO could well be a dP ferromagnetic 
semiconductor due to a subtle interaction between the 
spinal C2 ppTT antibonding state and the close-lying bot- 
tom of the host conduction band, as supported by hybrid 
functional calculations. Wc thus propose a mechanism — 
impurity pairing and molecular magnetism in resonance 
with the host band or not — to determine if a dP mag- 
netism could be present in C or N doped oxides. Ob- 
viously, this mechanism applies also to other open-shell 
impurities than those discussed in the present paper. We 
trust that the nature of the defect complexes giving rise 
to magnetism in MgO and ZnO will soon be identified by 
experiments. 
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